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ABSTRACT 

We present a simple theoretical framework for massive galaxies at high redshift, where the main 
assembly and star formation occurred, and report on the first cosmologies! simulations that reveal 
clumpy disks consistent with our analysis. The evolution is governed by the interplay between smooth 
and clumpy cold streams, disk instability, and bulge formation. Intense, relatively smooth streams 
maintain an unstable dense gas-rich disk. Instability with high turbulence and giant clumps, each a 
few percent of the disk mass, is self-regulated by gravitational interactions within the disk. The clumps 
migrate into a bulge in ,$10 dynamical times, or ,$0.5Gyr. The cosmological streams replenish the 
draining disk and prolong the clumpy phase to several Gigayears in a steady state, with comparable 
masses in disk, bulge and dark matter within the disk radius. The clumps form stars in dense 
subclumps following the overall accretion rate, ~ lOOJWoyr -1 , and each clump converts into stars 
in ~ 0.5 Gyr. While the clumps coalesce dissipativcly to a compact bulge, the star-forming disk is 
extended because the incoming streams keep the outer disk dense and susceptible to instability and 
because of angular momentum transport. Passive spheroid-dominated galaxies form when the streams 
are more clumpy: the external clumps merge into a massive bulge and stir up disk turbulence that 
stabilize the disk and suppress in situ clump and star formation. We predict a bimodality in galaxy 
type by z ~ 3, involving giant-clump star- forming disks and spheroid-dominated galaxies of suppressed 
star formation. After z ~ 1, the disks tend to be stabilized by the dominant stellar disks and bulges. 
Most of the high-z massive disks are likely to end up as today's early-type galaxies. 
Subject headings: galaxies: elliptical and lenticular — galaxies: evolution — galaxies: formation — 

galaxies: mergers — galaxies: kinematics and dynamics — galaxies: spiral — stars: 

formation 



1. INTRODUCTION 

The common picture of galaxy formation assumes that 
disk galaxies form by gas accretion and they transform 
into spheroidal stellar systems mostly by major merg- 
ers (Toomre & Toomre 1972; White & Rees 1978; Fall 
1979; Blumenthal et al. 1984; Mo et al. 1998). The disks 
are assumed to be the sites of quiescent star formation 
while the mergers are responsible for intense starbursts 
as well as subsequent quenching of star formation in the 
spheroidal merger remnant (e.g., Hopkins et al. 2007). 
While major mergers do occur, recent high-redshift ob- 
servations and theoretical developments indicate that 
this is not the major mode of galaxy formation and the 
star formation in them (see below). 

It appears that the most effective star formers in the 
Universe were galaxies of baryonic mass ~ 10 11 M Q in 
the redshift range z~ 1.5—3. The typical cases are repre- 
sented by UV-selected galaxies termed BX/BM (Adcl- 
berger et al. 2004) and rest-frame optically selected 
galaxies termed sBzK (Daddi et al. 2004), which we 
refer to collectively as massive "Star-Forming Galax- 
ies" (SFGs). Their mean comoving space density is 
n~2xl0~ 4 Mpc~ 3 (Tacconi et al. 2008), implying within 
the standard ACDM cosmology that they reside in dark- 
matter halos of masses M v < 3.5 x 10 12 M Q . The SFGs 
show star-formation rates of order 100M©yr _1 (Gen- 
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zel et al. 2006; Forster Schreiber et al. 2006; Elmegreen 
et al. 2007; Genzel et al. 2008; Stark et al. 2008), much 
higher than the few M Q yr _1 in today's Milky Way, while 
their baryonic masses and dynamical times are compara- 
ble to the Milky Way's. However, as argued below, their 
observed properties are incompatible with being ongo- 
ing mergers (Shapiro ct al. 2008) or merger remnants 
(Bournaud & Elmegreen 2009). The kinematics of many 
of the SFGs is consistent with a rotating disk, of circu- 
lar velocity V ~ 200 km s -1 and high velocity dispersion 
a ~ 50 km s~ 1 . The typical morphology is of a thick, gas- 
rich disk, extending to a radius i?d ~ 10 kpc. A unique 
feature of these galaxies, distinguishing them from low- 
rcdshift galaxies, is that the disk tends to be broken into 
several giant clumps, of ~ 1 kpc in size and up to a few 
times 10 9 M Q each, in which most of the star formation 
occurs (van den Bergh et al. 1996; Elmegreen et al. 2004; 
Elmegreen & Elmegreen 2005; Forster Schreiber et al. 
2006; Genzel et al. 2008). This is why they are sometimes 
referred to as "chain" or "clump cluster" galaxies. There 
is evidence that these clumps have formed internally in 
the disk (Bournaud et al. 2008; van Starkenburg et al. 
2008; Shapiro et al. 2008). In many cases there is also 
a significant central bulge, of a somewhat older stellar 
population than in the disk clumps (Genzel et al. 2008), 
but somewhat younger than in regular spirals (Elmegreen 
et al. 2008c). 

In contrast, an ongoing merger is not expected to show 
the kinematics and morphology of a rotating disk. The 
enhanced SFR in a merger is expected to occur either in 
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the two progenitor centers or in the compact coalescing 
core — not in the outer parts of an extended disk. Merger 
simulations show that a very wet merger with a special 
choice of orbital parameters may leave behind a non- 
negligible disk, but in most cases it is embedded in a more 
massive stellar spheroid remnant (Springel & Hcrnquist 
2005; Robertson et al. 2006; Robertson & Bullock 2008; 
Governato et al. 2009; Bournaud & Elmegreen 2009). 
These disks are expected to be smooth and not frag- 
mented into giant clumps (Robertson & Bullock 2008; 
Bournaud & Elmegreen 2009). Independently, the space 
number density of SFGs is higher by about a factor of 4 
than the expected density of merger-induced starbursts 
of comparable SFR (Elmegreen et al. 2007; Jogee 2008; 
Dekcl ct al. 2009). It should be noted that the most ex- 
treme star formers at high redshift are the dusty, bright 
Sub-Millimeter Galaxies (SMGs) (Chapman et al. 2004; 
Tacconi et al. 2008; Wall ct al. 2008), with SFRs of hun- 
dreds of M0yr -1 , but these are rarer cases, with space 
densities that are lower by an order of magnitude. While 
a significant fraction of the SMGs could indeed be star- 
bursts induced by major mergers, the much more com- 
mon SFGs pose the interesting question of how could 
such massive galaxies form most of their stars so rapidly 
at early times and through a process other than a major 
merger. A necessary condition is a steady intense gas 
supply. 

In parallel to the existence of SFGs, observations in- 
dicate that the other half of the galaxies of ~ 10 11 M Q 
at z ~ 2 are actually compact spheroids of radii ~ 1 kpc 
in which the SFR is suppressed to below ~ 10M Q yr _1 
(Kriek et al. 2006; van Dokkum et al. 2008). It seems 
that a division of galaxy type into a Blue Cloud and 
a Red Sequence, reminiscent of the robust bimodality 
known from lower redshifts (a review in Dekel & Birn- 
boim 2006), is already established by z ~ 2. A well devel- 
oped Red Sequence at z ~ 2 is not reproduced by current 
semi-analytic models of galaxy formation, and cannot be 
explained by the infrequent major mergers. The open 
questions are how do so many massive spheroids form so 
early, why are they so compact, and what is the mecha- 
nism responsible for the suppression of star formation in 
these galaxies. 

The gravitational fragmentation of gas-rich, thick tur- 
bulent disks into big clumps, and the subsequent migra- 
tion into a central bulge, have been proposed (van den 
Bergh et al. 1996; Elmegreen & Elmegreen 2005; Gen- 
zel et al. 2008; Bournaud et al. 2008) and successfully 
simulated for idealized disks in isolation (Noguchi 1999; 
Immeli et al. 2004b, a; Bournaud et al. 2007). Accord- 
ing to Toomre (1964), a disk becomes unstable once its 
surface density is high enough and the circular velocity 
and velocity dispersion are sufficiently low. The former 
is responsible for the local self-gravity and the latter tend 
to balance gravity against local collapse. Being gaseous 
helps maintaining a disk configuration because the gas 
cools while a stellar disk tends to acquire and maintain 
a high velocity dispersion and thicken with time. If the 
disk is marginally unstable, and if the velocity dispersion 
to rotation ratio is high, which is equivalent to having a 
high disk mass compared to the total mass within the 
disk radius, the typical clumps arc big. Being massive, 
the timescale for clump migration into the center, by 



clump-clump gravitational interactions and by dynam- 
ical friction, is short. The isolated simulations showed 
that the clumps coalesce into a central "classical" bulge 
(Elmegreen et al. 2008a) in several disk dynamical times, 
on the order of a few hundred Myr, leaving behind a sta- 
ble low-surface-density disk. 

However, we deduce from the appearance of giant 
clumps in a significant fraction of the massive galaxies at 
z ~ 2 (e.g., Elmegreen et al. 2007; Tacconi et al. 2008) 
that in each of these galaxies the clumpy disk configu- 
ration should last at least for a period comparable to 
the age of the Universe at that epoch, ~ 3 Gyr, rather 
than be a single short episode of migration. The key 
for making it a long-term phenomenon is, again, a con- 
tinuous, rapid supply of smooth cold gas that would re- 
plenish the disk as it is being drained. The incoming 
gas should keep the disk at the high density that ensures 
instability and new clump formation in the disk as the 
older clumps migrate inward. Indeed, such a continuous 
gas supply is very natural when the galaxy is consid- 
ered in a cosmological context at high redshift. The- 
oretical work (Birnboim & Dekel 2003; Dekel & Birn- 
boim 2006), and hydrodynamical cosmological simula- 
tions (Keres et al. 2005; Dekel & Birnboim 2006; Ocvirk 
et al. 2008; Dekel et al. 2009), have demonstrated that 
the galaxies in dark-matter halos of ~ 10 12 M© at z ~ 2 
arc typically Stream- Fed Galaxies, which arc built by a 
few cold narrow streams at rates ~ lOOM0yr _1 . Dekel 
et al. (2009) have found that, on average, two thirds of 
the stream mass are in rather smooth flows involving only 
clumps smaller than ~ 1O 1O M0, leading to mini-minor 
mergers of mass ratio below 1:10 relative to the existing 
galaxy. This has been shown to explain the abundance 
of high-rcdshift SFGs forming stars at high rates, three- 
quarters of which not through mergers. We show below 
that the rather smooth streams allow the continuous re- 
plenishment of a disk with high surface density without 
over-producing a velocity dispersion that could have sta- 
bilized the disk or completely destroyed it. 

On the other hand, we show that the clumpy com- 
ponent of the incoming streams does stir up the disk 
and can create a high velocity dispersion. Combined 
with growing a massive bulge, this may stabilize the disk 
against the formation of in situ giant clumps. If the 
high SFR is in the disk clumps, the elimination of these 
chimps is a quenching mechanism that suppresses the 
star formation preferentially in galaxies that are fed by 
clumpy streams and have a high bulge-to-disk ratio. 

We thus propose a scenario where the formation of 
galaxies is driven by cold streams, violent disk instabil- 
ity and the growth of a spheroid. The cosmological cold 
streams play a major role in this process. When their 
smooth component dominates, they build up the disks 
and make them fragment into giant clumps. These in 
situ clumps give rise to rapid star formation as well as to 
the formation of compact bulges. In the cases of clump- 
rich incoming streams, the external clumps help growing 
a massive spheroid through mergers and thus conspire 
with this spheroid to stabilize the disk and quench star 
formation. We argue that these are the two main modes 
of galaxy formation, dominating the buildup of massive 
galaxies at high redshift, where most of the mass has as- 
sembled into galaxies and most of the stars have formed. 



Massive Galaxies at High Redshift: Streams, Clumps & Spheroids 
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In this paper we present the theoretical framework for 
this scenario, via simple analytic estimates of character- 
istic timescales. In §2 we formulate the basic quantities 
relevant for giant-clump disk instability in terms of the 
key parameter 5, the ratio of disk to total mass within 
the disk radius. In §3 we address the self-regulation of 
the disk instability by the in situ clumps themselves. In 
§4 we estimate the timescales for clump migration and 
disk evacuation. In §5 we discuss the cosmological gas 
input rate through streams, smooth and clumpy. In §6 
we work out the evolution toward a steady state with in 
situ giant clumps when fed by rather smooth streams. In 
§7 we address the effect of clump-rich streams in growing 
a massive spheroid and joining forces with it toward sta- 
bilizing the disk and suppressing the SFR. In §8 we refer 
to the SFR within the disk clumps. In §9 we address 
the survivability of the in situ clumps, against rapid col- 
lapse and against disruption by stellar feedback. In §10 
we explain the spatial extent of the clumpy disks and 
the compactness of the spheroids. In §11 we address the 
instability of a combined disk if gas and stars and the dif- 
ference between the disks at high redshift and at low red- 
shift. In §12 we display preliminary maps from the first 
cosmological simulations that reveal disks in a steady- 
state giant-clump phase. These high-resolution simula- 
tions arc analyzed in a companion paper (Cevcrino ct al. 
2009). In §13 we summarize our analysis and discuss its 
implications. 

2. DISK INSTABILITY: GIANT CLUMPS 

2.1. Instability Criterion 

According to the standard Toomre instability analysis 
(Toomre 1964; Binney & Tremaine 2008, Chapter 6), a 
thin rotating gaseous disk becomes unstable to axisym- 
metric modes once the local gravity overcomes both dif- 
ferential rotation and pressure due to turbulence or ther- 
mal motions. This is expressed in terms of the stability 
parameter Q being smaller than a critical value of order 
unity, 



Q = ^< 



(1) 



Here ay is the radial velocity dispersion in the disk (or the 
gas sound speed if it is larger), £ is the surface density of 
the disk, and k is the epicyclic frequency. The latter is 
related to Q, the angular circular velocity at radius r, by 
k 2 = r dVL 2 jdr + 4ft 2 . The value of k ranges from f2 for 
Keplcrian orbits, through s/2Q for a flat rotation curve, 
in a uniform disk, to 2£1 for solid-body rotation. We 
adopt hereafter k = \/3Q, appropriate for high-redshift 
disks. Note that f2(r) is determined by the total gravi- 
tational force at r, which can be partly exerted by the 
disk itself and partly by a more spheroidal mass compo- 
nent. A value of Q below unity guarantees that there are 
scales of perturbations that are both Jeans unstable and 
rotation unstable, namely, they are large enough not to 
be stabilized by pressure (A > Aj = af/G£), and small 
enough not to be stabilized by the centrifugal force act- 
ing in the frame of the local perturbation due to the disk 
rotation (A < A rot = Att 2 GJ^/k 2 for a cold disk). 

For a thick disk, the analysis of radial modes is quali- 
tatively similar, as long as the perturbation length scale 
is smaller than the radius of the disk and larger than its 
thickness. It yields a slightly smaller value of the critical 



value, Q c ~ 0.68, for an isothermal thick disk (Goldreich 
& Lyndcn-Bcll 1965a). The pressure in the high-redshift 
thick disks, where ay ~ 50kms _1 , is clearly dominated 
by macroscopic motions, which we crudely refer to as 
"turbulence" . 

For a stellar disk, the instability criterion is similar to 
a gaseous disk, except that the factor n in the definition 
of Q is replaced by 3.36 — a negligible difference com- 
pared to the other uncertainties in the analysis. Thus, 
as long as the velocity dispersions of gas and stars are 
comparable, the instability analysis is valid for the com- 
bined system of gas and stars as a whole. We adopt this 
approximation for our simplified analysis of the high- 
redshift, gas-dominated disks. Observational evidence 
for this comes so far only from the thickness of the old 
stellar disk in "chain" galaxies as indicated by IR mea- 
surements, which is comparable to the thickness in the 
UV, associated with stars in formation and therefore gas 
(Elmegreen et al. 2009). In §11, we address deviations 
from this simple case, where we discuss star-dominated 
disks at lower redshifts. 

2.2. Clumps and Transient Features 

The characteristic wavelength of the unstable mode is 
(e.g., Binney & Tremaine 2008, Fig. 6.13) 



A,. 



2tt 2 G£ 



(2) 



This is ~ 0.5 A ro t and ~ 2Aj for Q ~ 1. The ring- like 
density perturbation is expected to break into lumps, 
and we take the pre-collapse radius of a typical lump to 
be R c ~ Ac/4. 

If the timescale for decay of pressure support in the 
growing perturbations is comparable to or shorter than 
the disk crossing time O -1 , the perturbations become 
gravitationally bound virialized clumps of radii ;$ R c . If 
the pressure decay time is longer, the perturbations are 
stretched by shear on a dynamical timescale and become 
elongated transient features on a length scale comparable 
to the disk scale (Goldreich & Lynden-Bcll 1965b). For 
example, Gammie (2001) simulated the case of a thin 
Keplerian gaseous disk and found that a cooling time 
as short as t coo \ < 3f2 _1 permits the fragmentation into 



bound clumps. He demonstrated that for t 



cool 



3fr 



the system reaches a gravo-turbulent steady state dom- 
inated by transient features in which the cooling is bal- 
anced by the dissipation of turbulence, by shocks or via 
a turbulent cascade to the viscous scale. In our case of 
a gas-dominated high-cr galactic disk, the gas radiative 
cooling time is much shorter than the dynamical time 1 . 
We thus expect a Q < Q c unstable disk to fragment into 
bound clumps encompassing a certain fraction a of the 
disk mass. Another fraction of the disk mass is expected 
to be in transient arm-like density perturbations, partly 
due to non-axisymmetric modes that could be unstable 
even for Q > 1. If the turbulence dissipation somehow 

1 The radiative cooling time is t coo i ~ 2.6 x 10 3 n~ 1 T4A~22 vr > 
where n is the gas density in atoms per cm 3 , here n ~ 10, T4 is the 
temperature in 10 4 K, and A_22(T) is the cooling rate in units of 
10 _22 Kcm 3 s _1 . At T4 ^ 1.8, A22 ~ 1, and at lower temperatures 
it drops sharply: at T4 = 1, A22 ~ 0.04, and below it is roughly 
A oc T 2 . Thus, for any T > 100K, the cooling time is much shorter 
than the dynamical time f2 _1 ~ 50Myr. 
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becomes slower than the dynamical time, e.g., because 
of fragmentation and star formation, the balance shifts 
from bound clumps to transient arms. In particular, the 
old stellar component, in which a does not decay, is ex- 
pected to be part of the transient features. 

Indeed, our simulations (§12) confirm the visual im- 
pression from the observed "chain" and "clump-cluster" 
galaxies that many of the disk clumps have well-defined 
round boundaries and are gravitationally bound physical 
entities. The simulations also show that the clumps are 
embedded in a perturbed disk with elongated transient 
features, and that a significant fraction of the clump mass 
is being exchanged with the surrounding arms and disk 
due to tidal effects. The notion that the clumps are long 
lived while the arms arc transients is supported by the 
finding that the high-rcdshift star-forming regions in spi- 
ral arms are younger than the stellar populations in the 
clumps of clump-cluster and chain galaxies (Elmegreen 
et al. 2009). 

2.3. Disk Fraction 

We adopt as our basic parameter the fraction of mass 
in the disk component within the characteristic radius of 
the disk Rd, 

r Md , . 

6 = ~~M ' (3) 

■'"tot 

where the total mass M to t within Rd includes the con- 
tributions of the disk and the spheroid of dark matter 
and stars. The maximum possible value of 5 is /3, the 
fraction of baryons including disk and bulge within the 
disk radius, 

6<P^^. (4) 

The ratio of disk to total baryonic mass is then 
Md/Mbar = 0~ X 5, so that a bulge-less disk is 5 = (3, 
and a disk- less bulge is <5 = 0. 
Expressing the circular velocity at Rd as 

(fii? d )W 2 c^, (5) 

and using the approximate relation 

we obtain a simple expression for Q: 



(6) 



(7) 



where we have assumed an isotropic three-dimensional 
velocity dispersion a = \[?>a r . The initial clump radius 
becomes 



6 



:6R d 



so the typical clump mass is 



36 S M *- 



(8) 



(9) 



It would be worthwhile evaluating the likely range of 
the disk fraction 5 and its maximum value f3 in the SFGs. 
We note up-front that the observed clumpy disks, with 
Q ~ 1, rule out naked disks, 5 — 1. This would have re- 
quired an extremely high dispersion, a ~ V, more com- 
patible with a spheroid. It would have also implied that 



the clumps are more massive than the observed clumps. 
Indeed, (3 is guaranteed to be less than unity due to the 
dark-matter contribution to the mass within the disk ra- 
dius. 

What is the typical value of the baryonic fraction (31 
While (3 ~ 0.5 in the Milky Way today, current estimates 
for SFGs at z ~ 2 range from 0.5 to 0.8 (Bournaud et al. 
2008; Forster Schreiber et al. 2009). One way to crudely 
estimate (3 from a theoretical perspective is by relating 
the disk radius to the halo virial radius R v via 



Rd ~ Ai? v 



(10) 



where A is roughly the halo spin parameter (Fall & Ef- 
stathiou 1980; Mo et al. 1998). The average value of A 
from cosmological TV-body simulations and tidal-torque 
theory is ~ 0.04 and rather independent of mass and 
time (e.g., Bullock et al. 2001a). On the other hand, 
the observed typical radii of Rd ~ 10 kpc for the high- 
rcdshift clumpy disks and the indicated virial radii of 
i? v ~ 100 kpc for the halos of these galaxies imply an ef- 
fective value of A ~ 0.1 for these galaxies (Bouche et al. 
2007; Genzel et al. 2008) (see a discussion of this exces- 
sive disk sizes in §10). Assuming further, very crudely, 
that the halo is an isothermal sphere of virial mass M v 
hosting a baryonic mass M om - of gas and stars, we can 
approximate M to t — XM V + M oar , and obtain 



/b + A 



(11) 



where f = M aar /M v is the baryonic fraction within 
the virial radius. The value of /b could be lower than 
the universal value of ~ 0.16 because of mass loss, e.g., 
due to supernova-driven winds (Dekel & Silk 1986). For 
A ~ fb ~ 0.1 we obtain ~ 0.5, and with A and /b in the 
ranges 0.04—0.1 and 0.05 — 0.15, respectively, the value of 
(3 ranges from 0.33 to 0.75. Replacing the crude isother- 
mal sphere with a more realistic NFW profile of a low 
concentration parameter C ~ 4, as appropriate for halos 
of the relevant masses at z ~ 2 (Bullock et al. 2001b), 
we obtain inside 0.1i? v a value of ~ 0.63. The obser- 
vational indications and this crude theoretical estimate 
lead us to adopt (3 = 0.6 as our fiducial value. 

Eq. (7) with Q ~ 1 implies that disk fractions in the 
range 5 ~ 0.3 — 0.6 are consistent with the observed range 
of velocity dispersions a r /V ~ 0.17 — 0.35 (Genzel et al. 
2008). For such values of 5 in cq. (9), the typical disk 
clumps are predicted to involve a few percent of the disk 
mass each, as observed (Elmegreen et al. 2007; Genzel 
et al. 2008). 

If each clump has contracted by a factor of two from 
its initial radius into virial equilibrium, the internal one- 
dimensional velocity dispersion within the clump, given 
by al ~ (l/3)GM c /(0.5-Rc), can be written using eq. (8) 
and (9) as 



V 



(12) 



Since the radial velocity dispersion of the clumps relative 
to each other, based on eq. (7), obeys a r /V ~ 6Q/y/3, 
we conclude that a c — a r when Q ~ 1. 
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3. SELF- REGULATED CLUMPY DISK 

For a disk of a given surface density in a given poten- 
tial well, namely given £ and f2, the instability provides a 
feedback loop that can drive an unstable disk into a self- 
regulated, marginally unstable state, where Q ~ Q c and 
the clumps are as massive as they could be. This means 
that the velocity dispersion ay is kept at the maximum 
possible value for which Q is still below the critical value 
for rapid instability. This self-regulation is achieved if the 
instability and fragmentation process itself is the gener- 
ator of velocity dispersion, and if it is capable of doing 
so on a timescale comparable to the timescale of turbu- 
lence decay. Then, if ay is temporarily low such that Q 
falls below Q c , the fragmentation becomes more efficient 
and it drives Q up toward Q c . If a r overshoots to high 
values such that Q > Q c , the disk becomes stable, the 
fragmentation process is suppressed, and a is allowed to 
drop such that Q settles to a value just below Q c . 

The turbulence of a gaseous disk dissipates by shocks 
and by a turbulent cascade to the viscous scale. With 
a turbulence energy of E ~ (3/ 2) A/a;?, and for a dis- 
sipation rate of E ~ Maf/R e ddy where i? dd y is the 
radius of the largest eddy, the dissipation timescale is 
idis ~ (3/2)i? c ddy/c r . If i?oddy — Rc, and we denote the 
disk dynamical crossing time by 



(13) 



we obtain from eq. (8) 



(14) 

For Q of order unity, the dissipation timescale is on the 
order of the disk dynamical time, with no explicit depen- 
dence on a or S. In order to maintain a self-regulated 
marginally unstable state, the stirring up of turbulence 
must occur on a comparable timescale. 

In an unstable disk with Q < 1, the density perturba- 
tions grow as e' w '*, where to obeys the dispersion relation 
(e.g., Binney & Tremaine 2008, eq. 6.66) 



-2irGY,\k\ + a 2 r k 2 



(15) 



One can see that the wave number of the fastest growing 
mode is k = 7rG£/a 2 , and the growth rate of this mode 
is given by uj 2 = k 2 — (7rGS / a r ) 2 . This growth rate, 
which vanishes for Q = 1, becomes comparable to the 
dynamical time for values of Q slightly below unity, e.g., 
\w\ = ft for Q = l/S/2 (assuming k = \/3il). The system 
should maintain a value of Q in this ballpark slightly 
below unity in order for the fragmentation to react in 
time to variations in the turbulence level. 

What is the mechanism that stirs up the disk and main- 
tains the required velocity dispersion? Some turbulence 
is generated by feedback from stars and supernovae (§9). 
Turbulence is also generated by dense gas clumps that 
flow in as part of the cold streams, but this process is not 
self-regulated and may stabilize the disk (see §7). How- 
ever, the highly perturbed fragmented disk is capable 
of self-regulating itself by the gravitational interactions 
within it, without the help of external energy sources, 
and thus maintain the disk in the marginally unstable 
state. One particular mechanism of this sort for generat- 
ing velocity dispersion in the disk is the gravitational en- 



counters between the bound disk clumps (see Wada et al. 
2002; Agertz et al. 2008; Tasker & Tan 2009). We demon- 
strate below that these clump encounters by themselves 
may be sufficient for self-regulating the disk instability. 

The timescale for the clump encounters to generate the 
velocity dispersion required for a given value of Q can be 
estimated as follows. First, given that the disk thickness 
is roughly (a z /V)Rd, with a vertical velocity dispersion 
g z ~ cr/y/3, the spatial number density of clumps is 



aM d /M c 
TTR 3 d a z /V ' 



(16) 



where a is the fraction of the disk mass that is in clumps 
at a given time. Second, the cross section for encoun- 
ters that involve an energy change of order a 2 can be 
estimated by 2 

a ~ ir{GM c /a 2 ) 2 . (17) 

Then, the characteristic time for the encounters to gen- 
erate a specific energy change ~ a 2 is estimated by 
i onc ~ (nccra)" 1 . Using eqs. (5), (7) and (9), we obtain 



t onc ~2.1a- 1 Q 4 t d 



(18) 



Recall that this expression is valid for Q of order unity, 
but note the strong dependence of t cnc on slight devia- 
tions of Q from unity, with no explicit dependence on S. 
The kinetic energy of the gaseous clumps could be dissi- 
pated in head-on collisions, but the cross section for such 
collisions is smaller by a geometrical factor of order a few 
than the cross section for effective gravitational encoun- 
ters. We therefore expect the net effect of clump encoun- 
ters to be of stirring up the dispersion on the timescale 
estimated in eq. (18). 

Self-regulation could be obtained by clump encoun- 
ters alone if t enc ~ idis, i.e., a~ 1 Q b ~ 0.67 when com- 
paring eqs. (18) and (14). This is valid, for example, 
when Q ~ 0.67 and a ~ 0.2. Recall that Q is ex- 
pected to be slightly below unity in the self-regulated 
unstable state both because Q c ~ 0.68 is the critical 
value for instability of a thick disk and because such a 
value of Q allows perturbation growth on a dynamical 
timescale. Independently, preliminary observations and 
simulations (Elmcgrccn & Elmegreen 2005; Elmegreen 
et al. 2007; Ceverino et al. 2009) indicate that the value 
of a can range from 10% to 40%. A value of a in this 
ballpark could indeed be expected if the unstable wave- 
length along the tangential direction is comparable to the 
radial wavelength A c , and if clumps form at the regions of 
positive interference of the tangential and radial waves. 
We thus adopt Q ~ 0.67 and a ~ 0.2 as our fiducial case. 

If the dissipation timescale somehow becomes signifi- 
cantly longer than the disk dynamical time, the balance 
between dispersion generation by clump encounters and 
turbulence decay is obtained over a timescale longer than 
id- Such a slowdown of the dissipation rate may result 

2 If the clump internal velocity dispersion cr; nt was smaller than 
the clump-clump velocity dispersion a, the cross section for the 
most effective individual encounter was actually nR 2 , which is 
larger by a factor (<r/cri nt ) 4 , but the energy change in each en- 
counter was smaller by a similar factor. Therefore, eq. (17) is 
indeed the effective cross section for a total change of ~ a 2 in a 
series of such encounters. 
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from the development of inhomogencitics and fragmen- 
tation of the disk gas, through the increase in size of the 
largest eddies and the possible decrease in cross section 
for eddy collisions. It may also happen when the stellar 
component of the disk becomes substantial, as idis is in- 
versely proportional to the gas fraction in the disk. In 
this case of slower dissipation the self-regulation in the 
marginally unstable state can be achieved by clump en- 
counters with a slightly larger value of Q and/or a lower 
value of a. 

The gaseous component between the clumps is ex- 
pected to participate in the instability and fragmentation 
process and thus to share the stirrcd-up clump veloc- 
ity dispersion. This is especially true for the large-scale 
transient features, which participate in the stirring up 
process as they are sheared and stirred up themselves. 
The inclusion of such additional stirring-up effects may 
be analogous to increasing a in the analysis of clump 
encounters above. 

4. CLUMP MIGRATION AND DISK EVACUATION 

The same clump interactions that generate velocity 
dispersion, and the dynamical friction exerted by the rest 
of the disk on the giant clumps, also affect the systematic 
rotation velocity of the clumps and make them migrate 
into the center. Both the smooth and clumpy compo- 
nents of the disk are participating in this process, so a 
should be replaced by unity for the purpose of estimating 
the migration time. Since the disk just outside a given 
clump has systematically lower angular velocity, it drains 
angular momentum from the clump. Over a timescale of 
at enc , it would systematically reduce the angular velocity 
of the clump by a and therefore cause it to migrate on 
a timescale a(V/cr)t cnc . However, the inner disk pushes 
the clump systematically outward on a similar timescale 
. The net effect acts on a timescale longer by (V/a) than 
that of each side of the disk alone (Goldrcich & Trcmaine 
1980; Ward 1997) The migration time is thus 



a 



V 



to 



2.1Q 2 5~ 2 t d 



(19) 



Note that the migration is rapid when Q is low and when 
S is high, namely when the disk is massive and unstable, 
but the migration timescale does not depend explicitly 
on a. 

Given that only a fraction a of the disk is in the giant 
clumps, the evacuation rate by clump migration of a mass 
comparable to the entire disk mass is 



(20) 



and the timescale for the evacuation of the entire disk 
mass is 



tr 



a t 



mig 



10.5 oT^ Q 2 S~ 2 t d , 



(21) 



where a 2 = a/0.2. With Q = 0.67 and a dominant 
disk of 6 — 0.35, say, we have t m i s — 7.6 t d - Then 
with a = 0.2, the evacuation time is £ cva c — 38 id- 
With t d ~ 50Myr, we obtain t mig ~ 380 Myr and 
i cvac ~ 1.9 Gyr. These estimates arc consistent with 
the findings in simulations of isolated disks (Bournaud 
et al. 2007). When 8 is slightly smaller, the timescales 



are somewhat longer, consistent with observational esti- 
mates (Elmegreen & Elmegreen 2005). 

The angular-momentum transfer and mass flow in the 
perturbed disk involve several other processes beyond the 
clump-clump interactions and dynamical friction. For 
example, the clumps induce transfer of disk angular mo- 
mentum outward, causing much of the disk to accrete 
inward at a rate that is comparable to the mass inflow 
rate directly associated with the clump migration itself 
(Sari & Goldreich 2004). This accretion rate exceeds the 
migration rate by d/h ^ 1, where d is the distance along 
the disk at which density waves shock and dissipate and 
h is the vertical scale height of the disk. 

The sheared transient features, which are also present 
in the unstable disk, exert torques which also cause 
angular- momentum transfer from the inside out. The as- 
sociated mass inflow rate, based on Shakura & Sunyaev 
(1973), is 

M sheaI ~ 3 s ' 2 G-^Q^a , (22) 

where a is the dimcnsionlcss angular-momentum flux 
density. According to Gammic (2001), it is given by 



a = [(9/4)7(7 -i)m dil 



(23) 



where 7 is the adiabatic index in the equation of state. 
For 7 = 5/3, appropriate for a turbulent medium, and 
with iltdis — 3 providing the maximum effect, we have 
a ~ 0.13. The timescale for evacuating the entire disk by 
the inflow associated with the transient features becomes 



tgheax — 7.5 Q 2 S 2 t d 



(24) 



With Q ~ 0.67, this is three times longer than the disk 
evacuation time by clump migration, eq. (21), but with 
Q closer to unity the two evacuation times become com- 
parable. 

We thus conclude that by considering only the clump 
migration, ignoring the other processes causing mass in- 
flow, we obtain an order-of-magnitude estimate for the 
disk evacuation rate that can serve as a lower limit to 
the actual evacuation rate. 

5. COSMOLOGICAL INPUT RATE BY STREAMS 

The key to understanding the high abundance of mas- 
sive disks with giant clumps at high redshift, as well as 
the high SFR in them, is the continuous vigorous input 
of gas into these disks via cold streams. The average rel- 
ative accretion rate into halos of mass M in the standard 
ACDM cosmology is given to a good approximation by 



^c0.47(l + z)P 5 M 1 V 5 Gyr" 1 



(25) 



where (1 + 2)3 = (l + -z)/3 and M 12 = M V /W 12 M Q . This 
has been derived by Neistein et al. (2006) as a useful fit to 
an analytic prediction based on the EPS approximation 
and it has been shown to fit well the halo growth rate 
measured in the Millennium cosmological N-body simu- 
lation (Neistein & Dekel 2008; Genel et al. 2008). Dekel 
et al. (2009) have shown using cosmological hydrodynam- 
ical simulations that at high redshift the same expression 
approximates the average relative baryonic input rate by 
cold streams into the disks of galaxies of a given bary- 
onic mass, Mbar/A^bar- The actual accretion rates into 
individual galaxies range from three times lower than 
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the average to twice the average. A fraction 7 of the 
stream mass is in clumps of mass larger than O.lMbar, 
with an average of 7 ~ 0.33 and an effective range from 
7 <C 1 to 7 as high as ~ 0.6 (we denote 7.33 = 7/0.33). 
The rest is in a smoother form consisting of mini-minor 
gas-dominated clumps and smooth gas in unknown pro- 
portions because of limited resolution in the simulations. 

We wish to relate this accretion rate to the disk dy- 
namical crossing time, which can be written as id = 
48 t4s Myr. A dynamical time of 48 Myr is obtained 
for a disk radius i?d = lOkpc and a rotation veloc- 
ity V = 200 km s , which are the typical values of 
~ 10 n M Q SFGs at z ~ 2. Then, eq. (25) defines the 
accretion time, 

t acc =^~44r- 1 (z,M)i d , (26) 

A/bar 

where the explicit redshift dependence and the weak 
mass dependence are in the factor 

r(z,M)=U 8 (l + z)l 25 M° 12 15 . (27) 

Note that r ~ 1 for the typical SFGs at z ~ 2. As long 
as the relevant SFG disks are observed to have similar 
radii and velocities at different redshifts, namely simi- 
lar dynamical times, eq. (27) is useful for describing the 
redshift dependence of t acc . 

Alternatively, it could be useful to use — \Ry , as in 
eq. (10), if A is rather constant in some redshift range, as 
deduced from cosmological simulations and tidal-torque 
theory (e.g., Bullock et al. 2001a). If the characteristic 
disk rotation velocity is comparable to the virial velocity, 
we have id — Ai v , where t v = R v /V v is the halo virial 
crossing time. The cosmological relation between halo 
virial radius and virial velocity 3 implies that r of eq. (27) 
can be alternatively expressed as 

t(z,M) = \. 1 (1 + z)° 3 75 M° 2 15 , (28) 

where A 1 = A/0.1. Note that if A is constant in time, 
as opposed to id being constant, the redshift dependence 
of i acc is weaker. The two expressions for r, eqs. (27) 
and (28), coincide at z ~ 2 for = A.i = 1. Based 
on the observed SFGs, it seems that the disk dynamical 
time is roughly constant at lower redshifts z < 2. On 
the other hand, a constant A may be more sensible from 
a theoretical point of view at z > 2 (see a discussion of 
disk sizes in §10 below). 

At z = 2, when the cosmological time is in — 3.25 Gyr, 
the virial crossing time is t v ~ 0.48 Gyr and i acc — 
2.1 Gyr compared to the disk dynamical time of t d — 
48 Myr. 



3 For the standard ACDM cosmology (fi m = 0.28, C A = 0.72, 
h = 0.7), in its Einstcin-deSitter regime, z > 1, the stan- 
dard virial radius and velocity are related to the virial mass as 

R v ~ 308kpc(l + ^)- 1 M 1 1 2 /3 and V v ~ H8kmB- 1 (l + «) 1 / 2 M^ 3 . 
Thus V200 — R100 (1 + z );s i where the quantities are in units of 
200kms~ 1 and 100 kpc, respectively. The virial crossing time is 
roughly a constant fraction of the Hubble time, t v ~ 0.15 in. In 
this regime the inverse of the expansion factor can be approximated 

by (1 + z) — 6.6tjj 2//3 , where the Hubble time tjj is in Gigayears. 



6. COSMOLOGICAL STEADY STATE 
6.1. Evolution Equation 

The high number density of massive SFGs that appear 
to be rotating disks broken into giant clumps indicates 
that they constitute a large fraction of the galaxies in 
halos of ~ 10 12 Mq at z ~ 2. This implies that the 
presence of giant clumps is not a short episode in the 
galaxy lifetime that disappears in one migration time, on 
the order of 400 Myr, but rather a long-term phenomenon 
that lasts for a Hubble time, on the order of 3 Gyr or 
more. The key for a long-lived disk with giant clumps is, 
again, the continuous, intensive gas supply into the disk 
by cold and rather smooth streams, valid at high redshift 
(Dekel et al. 2009). 

We wish to work out how 8 evolves in time, as the disk 
is being built and drained and as the bulge grows. The 
time derivative 8 can be derived from 



The baryon growth rate Afbax is the total baryonic ac- 
cretion rate from eq. (26), and the rate of change of disk 
mass, Md, is the difference between input to the disk and 
output from it, 

M d ~ (1 - 7) Af bar - M evac (8) , (30) 

with M evac ((5) from eqs. (19) and (20). For completeness, 
the associated growth rate of the bulge is 

M sph ~ 7 M bar + M evac (8). (31) 

The factor 7 allows for the possibility that only the 
smooth fraction (1 — 7) of the incoming streams adds 
to the disk mass, while the more massive clumps along 
the streams merge into the spheroid independent of the 
disk instability; that is, a standard merger. One assumes 
here that, on their way to join the bulge, those massive 
external clumps do not spend a significant amount of 
time as proper disk members (with a circular velocity 
V and a dispersion a) compared to the rest of the disk 
material. They therefore do not increase the disk surface 
density and thus do not speed up the migration and evac- 
uation of the in situ clumps. This assumption is justified 
because the massive external clumps are several times 
more massive than the internal disk clumps, and their 
effective masses are even larger because they are likely to 
carry with them dark matter, so their inward migration 
time by dynamical friction is much shorter than the disk 
evacuation time. From the simulation analyzed by Dekel 
et al. (2009), we learn that the average value of 7 ~ 0.33 
refers to external clumps more massive than 0.1Mb ar . 
This could perhaps serve as a crude estimate for the di- 
viding line between "massive" stream clumps that lead 
to major or minor mergers without joining the disk, and 
the smoother component involving "mini-minor" clumps 
that do join the disk and spend a non-negligible time 
there. We estimate further, based on the simulation, that 
about a third of the galaxies are fed by smooth streams 
of 7 -C 1, and another third have an excess of massive 
incoming clumps, 7 jt 0.5. 

As a simple example, we work out the case (3 = const., 
following the notion that the virial baryonic fraction /b 
and the effective spin parameter A that enter eq. (11) 
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are constant in time. The gas input rate is assumed to 
follow the cosmological mean for that mass and redshift. 
We also assume for simplicity that M v remains in the 
broad range 10 11 — 10 12 5 M Q , taking advantage of the 
very weak mass dependence of the accretion time and 
the dynamical time. Using the definitions cqs. (21) and 
(26) in eq. (30), we obtain for the rate of change of 8 



8 * p (1 - 7 - /T 1 *) t~ c \ - 6t- v \ c (6) 



(32) 



in which the ^-dependence is explicitly specified. We 
intend to solve this differential equation for 8(t). 

6.2. Steady State Attractor 
We first note that for a given accretion rate, there 
should be a steady-state configuration where 8 = 0; that 
is, 8 remains fixed at a critical value, S ss . This configura- 
tion serves as an "attractor" : if 8 > 8 SS , the efficient disk 
evacuation and spheroid growth make 8 decline toward 
8 SS , and if 8 < 8 SS , the efficient streaming into the disk 
forces 8 up back to 8 SS . In order to estimate the value of 
8 SS , we set 8 to zero in eq. (32), namely 



/?(l- 7 )<5- 1 -l 



'ovac 

The relevant time ratio, from eqs. (21) and (26), is 

^acc 



1 X2 



to 



4.18 a. 2 Q~ z t~ l 8 



(33) 



(34) 



so, the equation is a depressed cubic polynomial equation 
for 8 SS , 

8i + b8 ss -c = 0, 



b = 0.24a^ 1 Q 2 T , c= 6/3 (l- 7 ). 
The solution (e.g., by Cardano's method) is 

b 



3u 



C\l/3 

2 



46 3 



1/2' 



-i 1/3 



27c 2 J 



(35) 
(36) 

(37) 
(38) 



With the fiducial values Q = 0.67, a = 0.2 and (3 = 0.6, 
and with 7 in the range — 0.33, we obtain at z = 2, 



0.32- 0.25. 



(39) 



Simple solutions of eq. (35) in certain limits arc discussed 
in §A. 

The redshift dependence of 8 SS is rather weak. For 
example, with 7 = 0, we obtain 8 SS ~ 0.38 and 0.25 
at z = 9 and 1, respectively (where we have assumed 
A ~ 0.1 at z > 2 and t d ~ 48Myr at z < 2). The 
mass dependence is even weaker. Varying (3 in the range 
0.5 - 0.75 yields S ss ~ 0.28 - 0.35. Varying a in the 
range 0.1 — 0.4 gives S ss = 0.37 — 0.26. We see that when 
varying the parameters within their likely ranges about 
the fiducial values, the value of S ss at z ~ 2 remains quite 
stable about 0.3. Thus, the steady-state configuration is 
expected to be with a bulge mass comparable to the disk 
mass, 8 ~ j3/2, to within a factor of two. 
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Fig. 1. — Steady state and bimodality. The solid blue curves 
show the evolution of the disk-to-total mass ratio 8, according to 
eq. (32), for the fiducial case (Q = 0.67, a = 0.2, fa = 0.6, t d = 
0.015 1 and an average cosmological inflow rate). Shown are two 
extreme cases in terms of the presence of massive clumps in the 
incoming streams: very smooth streams (7 = 0) and very clumpy 
streams involving mergers (7 = 0.6). Eq. (32) is solved starting 
at z = 9 with either a dominant disk (5 = 0.5) or a dominant 
bulge (S = 0.1). Independent of the initial configuration, as long 
as the disk is unstable the system converges within less than a 
Hubble time to a quasi-steady state where 8 is slightly above 8 BB 
(long-dashed magenta curves). The critical value for stabilization 
by incoming stream clumps, <5 a t a b of eq. (43), is shown as a dashed 
red curve for f c = 0.33 and v = 1.2. It implies that streams that 
are smoother than average (7 < 0.33) produce a long phase of 
an unstable giant-clump disk with a comparable bulge that can 
last till after z ~ 1. Streams that are more clumpy than average 
(7 > 0.33) stabilize the disk and lead to a bulge-dominated system 
with a less massive, smoother disk. Therefore, the actual solution 
for 7 = 0.6 deviates from the lower curve shown. 

6.3. Convergence to Steady State 

With the average accretion rate evolving on a timescale 
that is not much longer than the disk evacuation 
timescale, we need to verify that the system indeed ap- 
proaches a quasi-steady state, and find out how far it 
might be from <5 SS where 8 = 0. Figure 1 shows numer- 
ical solutions for 8(t) of eq. (32) with eq. (34), for the 
fiducial choice of parameters and id = 0.015 1, eq. (28). 
The initial conditions are set, as an example, at the reion- 
ization epoch, z- m ~ 9 or t- m ~ 0.54 Gyr. We choose the 
starting configuration to be alternatively disk dominated, 
<5 in = 0.5, or bulge dominated, 8 iu -C 1. Independent of 
the initial conditions, the system settles within a Hub- 
ble time to a value of 8 near 8 SS , which is also shown. 
This can be understood via the approximate solutions to 
eq. (32) discussed in §A. The actual steady-state value 
of 8 is slightly higher than i5 ss because the latter declines 
with time on a timescale that is somewhat larger than 
the timescale for changes in 8, driven in this range by 

tcvac • 

We conclude that throughout the broad redshift range 
of interest the system tends to be in a steady state in 
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which the bulge and disk masses are comparable, to 
within a factor of two. The disk fraction in this steady 
state is larger in cases where the fraction of massive 
clumps in the incoming streams is smaller. Note that 
the evolution curves shown in Fig. 1 are derived, so far, 
under the assumption that the disk is always unstable 
to the formation of giant clumps. We will see next that 
this assumption may be valid only for streams that are 
smoother than average, namely a low 7 and then a high 
6. An excess of incoming steam clumps may stir up the 
disk to stabilization with a low value of 6. 

The observed high-redshift SFGs indeed seem to indi- 
cate convergence into a steady state with a bulge of mass 
comparable to the disk mass, 6 ~ 0.3, as predicted here. 
Genzel ct al. (2008, and private communication) have 
measured for six SFGs the "concentration" C, defined as 
the ratio of total dynamical masses within the radii of 
0.4" and 1.2", which roughly correspond to (l/3)i?d and 
i?d, respectively. Allowing for 20% of the mass within 
the inner radius to come from dark matter and disk (ap- 
propriate for an NFW halo profile and a uniform disk 
profile), the observed value of 0.8C could be interpreted 
as the bulge-to-total mass ratio within the disk radius, 
0.8C ~ 13-8. When they plot in Fig. 5 the value of 
C against the time t* of onset of star formation in that 
galaxy 4 , they find that C is rising monotonically with 
t* for t* < 0.5 Gyr until it settles to a rather constant 
value near C ~ 0.38 for t* ~ 0.5 — 2 Gyr. This is con- 
sistent with the quick convergence to a steady state as 
predicted in Fig. 1, and with the steady-state configu- 
ration of 5 ~ 0.3 (for (3 = 0.6 and 7 < 0.33), which 
corresponds to a bulge-to-total ratio 0.8C ~ 0.3. 

7. STABILIZATION BY EXTERNAL STREAM CLUMPS 

Beyond the self-regulated stirring up of the disk by the 
in situ disk clumps themselves (§3), the incoming streams 
provide an alternative source for generating velocity dis- 
persion in the disk. Unlike the stirring by the internal 
products of the disk instability, which is naturally self- 
regulated, the external source could in principle push a 
up to a level that stabilizes the disk. 

When a smooth stream of gas density p s hits a disk 
of much higher gas density pd, conservation of momen- 
tum implies that the kinetic energy added to the disk 
is roughly a fraction p s / Pd of the total kinetic energy of 
the stream, while most of the incoming energy turns into 
heat that is radiated away. This implies that the smooth 
component of the streams, which is less dense than the 
disk by more than an order of magnitude, is not effective 
in stirring up the disk. 

However, the streams contain a clumpy component 
with gas densities comparable to the disk density. In the 
impact of such an external dense clump with the disk, 
about half the kinetic energy is expected to turn into 
turbulence in the disk. The timescale for stirring up the 
disk to a given value of a by external gas clumps that 
stream in with a velocity Vj n can thus be estimated by 

M d a 2 , . 

4 t„ is the time to form the current stellar mass at a constant 
SFR that equals the current SFR — a model that fits the data well 
(Forster Schreiber et al. 2009) 



The assumption here is that incoming clumps that con- 
stitute a fraction 7 C of the stream mass share most of 
their kinetic energy and angular momentum with the 
disk. This "f c is not necessarily the same as the 7 used in 
§6 to compute the evolution of an unstable system into 
a steady state. It may include the same massive clumps, 
under the assumption that they stir up the disk before 
they end up merging with the bulge, as well as mini- 
minor clumps (e.g., less massive than 0.1Mb ar ), which do 
not merge quickly with the bulge but are dense enough 
for stirring up the disk. We defer a more accurate esti- 
mate of the fraction of dense clumps in the streams to 
a future work based on simulations with higher resolu- 
tion. Meanwhile, for the purpose of the current simplified 
analysis, we adopt j c = 0.33, the average value deduced 
from the MareNostrum simulation used in Dckcl et al. 
(2009) 5 . We mark v = Vi n /V v , and read from the same 
simulation that on average v 2 ~ 2. 

We see from eqs. (19) and (21) that a 2 /V 2 ~ t cnc /t cv&c , 
so with tg.ee from eq. (26), eq. (40) becomes 



i stab ~2 7c - 1 U - 2 /?- 1 



St 



(41) 
(42) 



The first equality tells us that when the external clump 
accretion rate is sufficiently higher than the disk evac- 
uation rate, the stirring by external stream clumps is 
more effective than the stirring by internal disk clumps, 
which operates on a timescale ~ t cnc - This generation 
of turbulence by an external source is not self-regulated 
and could stabilize the disk if it operates on a timescale 
that allows it to balance the dissipation of turbulence, 
namely if i sta b < ^dis- Using eqs. (14) and (42), this 
condition translates to stabilization by external stream 
clumps once 



5<4tab^0.25( 7cU 2 /3r) 1 / 3 i 



T 1 • (43) 

Naturally, the disk tends to be stabilized by the streams 
when the streams contain more dense clumps and the 
accretion rate is high, and especially when d is low due 
a significant contribution from the spheroid. The latter 
is because a low <5 requires a low critical a for Q ~ Q Cl 
which makes it easier for a given stream to drive a above 
the critical value. 

Note that the mass dependence of <5 s t a b, entering 
through t in eq. (43), is very weak. For a halo mass 
smaller by an order of magnitude, <5 s tab is smaller by 
~ 10%. The redshift dependence entering through t is 
somewhat more pronounced, but it is still rather weak; it 
depends on whether id or A are assumed to be constant 
in time, eqs. (27) or (28), namely <5 sta b oc (1 + z) - 75 or 
(1 + z) 25 , respectively. 

Figure 1 shows <5 s tab for the fiducial case 7 C = 0.33 
and v = 1.2. We learn that for average dumpiness, 
7 ~ 7 C ~ 0.33, the value of 5 s tab happens to be in the 
ballpark of the steady-state value of S. The value of 
7 ~ 0.33 thus marks the transition between two very dif- 
ferent evolution tracks. In the cases of streams smoother 

5 It may be interesting to note that Khochfar & Silk (2009) esti- 
mate from a semi-analytic model that ~ 18% of the stream kinetic 
energy is required for driving a velocity dispersion as observed, con- 
sistent with the 7c /2 factor predicted by our theoretical argument, 
with our average value 7 C ~ 0.33. 



10 



Dckcl, Sari and Ccvcrino 



than average, 7 < 0.33, the stirring by external stream 
clumps is unimportant, as 5 > <5 s tab at all times. This is 
especially true for the case 7 « 1 shown in the figure. 
In this case, and as long as the accretion rate follows the 
cosmological average, the disk is in a long-term giant- 
clump phase near steady state with 6 slightly above S ss . 
For 7 « 1, the value of 6 ranges from 5 ~ 0.4 at very 
high z to 8 ~ 0.3 at z ~ 1, namely the bulge contains 
between a third and a half of the baryonic mass. This is 
compatible with the bulges detected in typical SFGs at 
high rcdshift, as mentioned at the end of §6. 

In the cases where the streams are more clumpy than 
average, 7 > 0.33, the stabilization by external stream 
clumps becomes important, as 6 < i5 s tab at all times. The 
disk does not develop a giant-clump phase (and if it some- 
how did at very high redshift, it was a short episode). 
In such a case, the bulge grows predominantly by the 
merging incoming massive clumps. The system does not 
evolve according to the evolution track shown in Fig. 1 
for 7 = 0.6. Instead, it settles to either S ~ <5 sta b or 
6 = (3(1 — 7), whichever is smaller. In this case the bulge 
is a high fraction of the baryonic mass, on the order of 
7- 

Our analysis thus predicts a robust bimodality in the 
galaxy population already at high rcdshift, driven by the 
spread in the degree of dumpiness in the streams feed- 
ing the galaxies. A small deviation of 7 from the av- 
erage value ~ 0.33 triggers a qualitative difference in 
the evolution pattern, even if the overall mass growth 
rate follows the cosmological mean. This bimodality is 
consistent with the observed high abundance of clumpy 
star-forming, extended disks at z ~ 2, alongside with 
the realization that a significant fraction of the massive 
galaxies are dominated by compact spheroids and show 
low SFR already at these redshifts (Kriek et al. 2006; 
van Dokkum et al. 2008). We predict that many of these 
high-rcdshift spheroids arc surrounded by smooth, sta- 
ble disks containing about a third of the baryonic mass, 
which continue to evolve secularly and form stars qui- 
escently (§11). Such secondary disks are indeed appar- 
ent in some of the HST NIC2 images of compact passive 
spheroidal galaxies at z ~ 2 — 3 (van Dokkum et al. 2008, 
Fig. 1). This configuration, of a dominant compact stel- 
lar spheroid and a smooth secondary disk, is not very 
different from the predicted remnants of certain binary 
wet major mergers (Springcl & Hcrnquist 2005; Robert- 
son & Bullock 2008; Governato et al. 2009), not to be 
confused with the extended, high-redshift disks of giant 
clumps forming stars that are not dominated by the cen- 
tral spheroid (Bournaud & Elmegrecn 2009). 

The evolution described in Fig. 1 assumes continuous 
gas input at the average cosmological rate and a constant 
level of dumpiness in the streams. Time variations in the 
stream properties can shift a galaxy from an unstable to 
a stable mode and vice versa. An example is illustrated 
in the simulated evolution of a Milky- Way-type galaxy 
termed "Via Lactea" ?, Fig. 5, where the overall growth 
rate drops to below-average levels between z ~ 2.5 and 
1.8. As a result, the bulge-to-disk ratio grows and i5 drops 
below 5 s tab, bringing the unstable giant-clump phase to 
an end. When the accretion resumes at z ~ 1.8, it takes 
the system a Hubble time, about 4 Gyr, to rebuild a disk 
that may bring it back to the unstable regime with 5 > 
<5stab after z ~ 1. By that time, the disk is already in the 



star-dominated regime and may never resume instability 
(§11)- 

8. STAR FORMATION 

The observed SFGs of M bar ~ 10 n Af Q at z ~ 2 show 
that the giant disk clumps are the sites of intense star 
formation. This indicates that the formation of clumps 
is a necessary condition for the main mode of star forma- 
tion in high-redshift galaxies. While the physics of this 
mode of star formation is beyond the scope of the current 
paper, we discuss here certain implications of our analy- 
sis of gas streams and disk instability, combined with the 
observed overall SFR, on the process of star formation 
in the clumps. 

The observed overall SFR of ~ 100M Q yr _1 across the 
disk (Elmegreen et al. 2007; Genzel et al. 2008) is only 
slightly smaller than the theoretically predicted gas in- 
put rate via cold steams into the central galaxies of dark- 
matter halos that are comparable in abundance to these 
SFGs (Dekel et al. 2009). This leads to the approxi- 
mation that the overall SFR in the disk is a significant 
constant fraction of the gas accretion rate from eq. (26): 
M*d — s Mb ar with s ~ 0.5. The stars formed earlier 
in the incoming external clumps are not part of the star 
formation in the in situ disk clumps. 

As a simple model for the SFR in the disk clumps, we 
assume that the clump mass M c = M g + M* remains 
constant as it turns gas into stars, and that the SFR in 
the clump is M* c = M g /t* with a constant t*. Then 

M* c (t) = (M c /i*) exp(— t/ttt). This means that most of 
the stars form during a period ~ t* with a SFR 011 the 
order of M c /t*. The fact that a significant fraction of 
the accreted gas is observed to form stars in the clumps 
while they are still in the disk implies that t* is not much 
longer than the clump migration time t mlg . On the other 
hand, the observations indicate that the number of star- 
forming big clumps, ~ 5, is comparable to the number of 
clumps seen in simulations (§12) and consistent with our 
fiducial values of M c and a from the instability analysis. 
If the clumps survive for a migration time (but see §9.2), 
this implies that t* is not much smaller than t m i s . We 
conclude that t* ~ i m ig- 

This is consistent with the cosmological steady state 
(§6.2) because the requirement that the SFR be compa- 
rable to the accretion rate is equivalent to the approxi- 
mate condition for steady state, i C vac ~ iacc- Indeed, if 
we write t* = aM^/ iW*<j, we obtain 

U-ap-h^dt^ (44) 

~8.8a. 2/ 9 _1 s _1 r _x <yid, (45) 

which is comparable to i m ; g from eq. (19) for our fiducial 
unstable disk with 5 = 0.3 at z ~ 2, where £» ~ 8.8 id ~ 
420 Myr. 

Based on the above simple model, the gas fraction 
in the clump when it coalesces with the bulge is / g = 
exp(— £ m ig/£«) — e _1 . This implies that the disk clumps 
merge into the central spheroid while they are still rela- 
tively gas rich. The process is thus similar to a wet ma- 
jor merger, where the gas shocks, heats up, cools, and 
condenses to the center, leading to a compact spheroid 
significantly smaller than the disk size (e.g., Dekel & Cox 
2006; Covington et al. 2008). This may explain the com- 
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pactness of observed spheroids at z ~ 2 (van Dokkum 
et al. 2008). In the late phases of spheroid growth, when 
6 is smaller, t m i g /t* is larger, so the mergers of the disk 
clumps into the spheroid are less dissipativc, making the 
late-formed spheroids less compact. 

The SFR efficiency in the disk clumps, based on 
eq. (45) and with t c ~ 0.5id, is 

' M c /t c U V ; 

~ 0.06 o^stJ -1 . (47) 

For the fiducial case, with 5 ~ 0.3 at z ~ 2, we obtain 
?7 ~ 0.06. Similar high efficiencies were estimated for 
z ~ 2 SFGs based on the pressure in the star-forming 
clumps (Lehnert et al. 2009). However, such values of r] 
are higher than the efficiencies in the most efficient star- 
forming clouds at low redshift. The common wisdom, 
based on the empirical Kennicutt-Schmidt law (Kenni- 
cutt 1989) and on star-formation theory (Krumholz et al. 
2009), is that r\ is of order 1% in all star-forming environ- 
ments, independent of the actual density. This indicates 
that the giant disk clumps may not form stars uniformly 
over their whole volume. Instead, stars form with a stan- 
dard efficiency ~ 1% in subclumps where the density is 
two orders of magnitude higher so the local dynamical 
time is an order of magnitude shorter. This is an observ- 
able prediction, to be addressed when the star-forming 
clumps are resolved. 

9. CLUMP SURVIVAL AND DISRUPTION 

The appearance of giant gas clumps forming stars in 
the disks of many SFGs indicates that they survive for 
at least several dynamical times. The same is evident 
in simulations (see §12). One wishes to understand the 
mechanisms that prevent the clumps from collapsing into 
themselves on a free-fall time and help them survive the 
various processes that work to disrupt them. In par- 
ticular, once the disk clumps form stars at high rates, 
they could be affected by stellar feedback, in the form 
of supernovae-driven winds or radiative feedback from 
massive stars, which we briefly comment on below. A 
preliminary discussion of the possible gravitational ori- 
gin for the turbulence pressure supporting the clumps is 
provided in §B. 

9.1. Supernova Feedback 

Could supernova feedback provide the pressure sup- 
port against free-fall collapse? If so, is there a risk that 
supernova-driven winds may remove the remaining gas 
and disrupt the clumps before they complete their mi- 
gration? Dekel & Silk (1986) evaluated the maximum 
energy fed into the interstellar gas by supernovae, tak- 
ing into account the radiative loses, £?sn — M*t T&< \, 
where M* refers to the SFR in each clump, e is the en- 
ergy released by a typical supernova (~ 10 51 erg), and v 
is the number of supernovae per unit mass of forming 
stars (which for a typical IMF is v ~ 1 per 50M©). The 
characteristic time i ra d marks the end of the "adiabatic" 
phase and the onset of the "radiative" phase of a typi- 
cal supernova remnant, by which it has radiated away a 
significant fraction of its energy and became ineffective. 
They found that it takes a similar time for the expand- 



ing supernova shells to reach a significant mutual over- 
lap, which allows an even distribution of the supernova 
energy (minus the radiative loses) over most of the gas. 
A necessary condition for the clump to be significantly 
affected by supernova feedback is that -Esn exceeds the 
binding energy of the clump, ~ (l/2)M c er 2 . This occurs 
if the potential well associated with the clump is shallow 
enough, with a three-dimensional velocity dispersion 

a*<V$ s = 2evr l ^. (48) 

Ed 

It turns out that in the relevant temperature range the 
cooling rate scales approximately as A oc T , which im- 
plies that the ratio of timescales is roughly t ra d/id ~ 0.01 
independent of gas density and other clump parameters. 
Dekel & Silk (1986) then assumed maximum efficiency in 
a burst of star formation, rj ~ 1, and obtained a critical 
velocity of Vsn — 100 kms -1 as an upper limit for the 
virial velocity of a system in which one may expect sub- 
stantial gas heating or removal by supernova feedback. 

With the estimated SFR efficiency of 77 ~ 0.1 in each 
clump as a whole, the critical velocity is reduced to 
Vsn ~ 30 kms -1 . This implies that supernova feedback 
cannot be effective in the giant clumps during the steady- 
state instability phase of a V ~ 200 kms -1 disk, where 
8 ~ 0.3 and therefore a ~ 40 — 60 kms -1 by eq. (7). 
The clumps may become more vulnerable to supernova 
feedback, a < 30 kms -1 , once 5 < 0.2, namely when 
the bulge becomes twice as massive as the disk (see also 
simulations by Tasker & Bryan 2008). However, when 
the clump fragments, the supernova-driven winds may 
escape harmlessly via low density "chimneys" through 
a porous interstellar medium (e.g., Ceverino & Klypin 
2009), and thus not lead to disruption even in clumps of 
a well below 30 kms -1 . 

9.2. Radiation Pressure 

While supernova feedback may be unimportant in the 
giant clumps, radiative feedback from O stars may be 
more effective. Murray et al. (2009) estimate that the 
clumps should be disrupted by radiative pressure once 
about 20% of their gas has turned into stars, expelling the 
remaining 80% of the gas back to the disk. If true, we can 
deduce that the system would maintain a steady state 
similar to the one described in §6, with the disk evacua- 
tion by migration replaced by an "evacuation" into dense 
star clusters that do not participate in the disk instabil- 
ity any longer. The two steady states are similar because 
the timescales for clump migration and for star formation 
in a clump mass are comparable (§8). The expelled gas, 
combined with the gas streaming in from the outside, 
would then help keeping the disk gaseous and unstable. 

However, if each clump is disrupted this way, its mass 
becomes smaller by a factor of ~ 5, so the bulge buildup 
by migration becomes slower by a factor of ~ 5 2 . This 
is because the migration time is inversely proportional 
to the clump mass, and the mass added to the bulge by 
each coalescing clump is proportional to the clump mass. 
This implies that most spheroids have to be made by 
mergers, in possible conflict with the cosmological merger 
rate and with most semi-analytic simulations of galaxy 
formation (also discussed in §13). On the other hand, 
this scenario predicts the common existence of massive 
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thick stellar disks, constructed from the slowly migrating 
star clusters, which may be associated with today's thick 
disks and SO galaxies. 

The stellar populations in the high-rcdshift clumps can 
provide an observational test that could clarify the sur- 
vivability of the clumps. With r] ~ 0.06, one fifth of the 
clump mass becomes stars in about three clump dynam- 
ical times, namely ~ 75 Myr. If radiative feedback sup- 
presses star formation after that time, the spread of stel- 
lar ages in each clump should not exceed ~ 50— 100 Myr. 
The slow migration implies that the single-age popula- 
tion in the clumps could be as old as ~ 2 Gyr. In con- 
trast, if the gas clumps survive for a migration timescale, 
the age spread could be ~ 0.5 Gyr, and giant clumps 
should not show populations older than ~ 1 Gyr. In this 
case, we predict that the typical giant clumps should be 
gas rich and forming stars at a high rate while they also 
contain a stellar population that has been formed in a 
similar rate during the preceding few Gigayears. There 
is observational evidence in favor of the latter (Forster 
Schreiber et al. 2009) Furthermore, if the clumps survive 
for a migration timescale, one expects most of the clumps 
at large disk radii to contain stellar populations younger 
than ~ 0.5 Gyr, while the clumps at smaller radii may 
also contain older populations. In contrast, if the clumps 
have lost their gas reservoir in less than ~ 100 Myr, they 
must have formed all their stars near the radius where 
they were formed and where they orbit for more than 
1 Gyr, thus showing no obvious age gradient with ra- 
dius. Finally, the massive radiation-driven outflows from 
clumps should be observable as systematic blueshifts at 
the clump locations. A preliminary search for such out- 
flows yielded a null result (K. Shapiro, private commu- 
nication). 

10. GALAXY SIZE 

As mentioned in §1, the observed SFG disks of Md ~ 
10 11 Mq at z ~ 2 are surprisingly extended, with radii 
~ lOkpc, comparable to the Milky Way today. These 
disk radii indicate A ~ 0.1 in eq. (10) whereas the average 
value is expected to be smaller by a factor of 2 — 3 (Bul- 
lock et al. 2001a). We note that the proposed scenario 
of stream-driven unstable disks provides several possible 
explanations for the extended disk sizes, as well as for 
the compact spheroids. 

The mass inflow in the disk into its center, by clump 
migration due to dynamical friction and encounters or 
by torques involving the transient shear and tidal fea- 
tures (§4), is associated with angular momentum trans- 
fer into the outer disk, which tends to stretch it further 
out. If the baryons in the bulge have lost all their angu- 
lar momentum to the disk, by the time the bulge mass is 
comparable to the disk mass, the angular momentum per 
unit mass in the outer disk should be about twice the av- 
erage, so the disk is expected to be twice as extended as 
it would have otherwise been. The timescale for reaching 
this state, eqs. (21) and (24), is on the order of one or two 
orbital times, namely a few hundred Megayears, and the 
configuration of comparable disk and bulge, 5 ~ 0.3, is 
indeed the typical configuration expected in the clumpy 
steady state (§6). 

Furthermore, in massive high-z disks that are fed by 
smooth streams, the generation of clumps and therefore 
the regions of intense star formation are expected to ex- 



tend to the outer disk. First, the streams that determine 
the angular momentum of the disk tend to be coherent 
and to join the disk with an impact parameter compara- 
ble to the disk radius. They therefore generate high gas 
surface density £ at large radii, as opposed to the expo- 
nential profiles of today's disks. The relative smoothness 
and low density of the high-z streams make them ineffi- 
cient drivers of turbulence (§7), so they are not expected 
to push a in the outer disk to high values. Finally, with 
a near flat rotation curve driven by the dark halo and 
the rather uniform disk outside the bulge, the angular 
velocity declines with radius roughly as f2 oc r -1 . These 
three factors allow Q oc crfi/E to become low and ad- 
mit instability at large radii. This is seen in simulated 
galaxies (Ceverino et al. 2009) and is consistent with the 
appearance of a broad, large-radius ring of giant clumps 
in some of the SFGs (Gcnzcl et al. 2008). 

The situation is very different in low-rcdshift disks. 
Unlike the halos at z ~ 2 that are fed by intense, co- 
herent, narrow streams, a 10 12 M Q halo at low redshift is 
fed from all directions by rather slow and erratic accre- 
tion (Keres et al. 2005; Dekel & Birnboim 2006; Ocvirk 
et al. 2008; Dekel et al. 2009). The exponential profile 
tends to provide a sufficiently high £ only at small radii, 
where the halo-driven rotation curve is closer to a solid- 
body rotation with a constant fl, so the instability is not 
preferred at large radii. 

Another possibly relevant mechanism is supernova 
feedback from early star formation in the central regions 
of galaxies. It is expected to be more effective at high 
redshift, when both the gas density and the SFR were 
higher. Supernova-driven winds may preferentially re- 
move low angular momentum gas that has settled at ear- 
lier times in smaller radii, and thus leave behind a more 
extended system with higher angular momentum (Mailer 
& Dekel 2002; Governato et al. 2007; Scannapieco et al. 
2008). 

Observations by Law et al. (2009) indicate that among 
the high-rcdshift galaxies with baryonic masses of a 
fewx 10 10 M Q , slightly below the range of SFGs addressed 
in our paper, about a third seem to be supported by 
velocity dispersion rather than rotation, with a/V ~ 
1/2 — 1. These galaxies tend to be gas rich (up to 90%) 
and not overly extended, with typical radii R ~ 1 — 2 kpc. 
Such a galaxy could be formed if fed by cold flows of low 
impact parameters that carry little angular momentum. 
Alternatively, if these flows are overly clumpy, they could 
be effective in stirring up an excessively high dispersion 
in the gaseous galaxy. The absence of a substantial, un- 
stable disk component in such a galaxy during its history 
may explain the relatively low past SFR and thus the ob- 
served high gas fraction. 

11. GAS & STELLAR DISKS AT HIGH & LOW REDSHIFTS 

The high-redshift disks are observed to be gas rich, at 
the level of 50% and perhaps more (Bouche et al. 2007; 
Daddi & et al., 2008). In contrast, today's massive disks 
are dominated by stars, with gas fractions of ~ 20% or 
less. Whereas both the gas and stellar components of 
the disk participate in the disk instability through their 
contributions to the local self-gravity, the gas dissipa- 
tion introduces a long-term difference between the two 
components. While the decay of gas turbulence and gas 
cooling allow the velocity dispersion or sound speed to 
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decline even when the system is self-regulated and in a 
cosmological steady state, the stars tend to maintain the 
velocity dispersion that they have acquired during their 
history. The young stars may still have a a that is not 
much larger than that of the gas, but the older stars are 
likely to have a higher a. Therefore, the young stars 
actively participate in driving the perturbation growth 
and they may even follow the gas into the bound clumps, 
while the older stars tend to form transient perturbations 
and eventually join the stable component that just adds 
to the external potential well and thus helps stabilizing 
the disk. 

The cosmological gas accretion rate plays an important 
role in determining the gas fraction in the disk. If stars 
form only in massive gas clumps, the disks are expected 
to remain gas rich as long as the gaseous input into the 
disk manages to replenish the disk mass on a timescale 
shorter than the time for the entire disk to turn into 
clumps, which is comparable to i C vac, namely 



1 < 



(1 - 7)^1 



ar 



0.24/3(1 -7) a 



"a 1 ' 



(49) 



The second equality is based on eqs. (21) and (26). For 
the unstable disk in steady state (§6.2), with 7 < 0.33, 
this ratio is ^ 2 — 3 at z ~ 2, and it varies slowly with red- 
shift. This implies that the disk should remain gas rich 
as long as the input is dominated by gas, which is more 
likely at high redshift. This indicates that the analysis 
performed in the current paper so far is a useful approxi- 
mation for the high-2 SFGs. At z < 1.5, the penetration 
of cold streams into the centers of halos more massive 
than 10 12 M Q becomes limited (Dekel & Birnboim 2006; 
Ocvirk et al. 2008; ?), so the input rate of cold gas to 
the disks is lower than implied by eq. (25). At these late 
epochs, the fraction of stars in the accreting matter be- 
comes higher. The disks gradually become dominated by 
the stars, in clumps and in transient features or spread 
out in the disk by shear and tidal stripping. 

A more accurate analysis of disk instability, espe- 
cially at low redshift, should therefore deal with multi- 
component disks. The axisymmetric instability of a two- 
component disk has been studied by Jog & Solomon 
(1984) and Rafikov (2001). Denoting the velocity disper- 
sions of stars and gas cr* and a g , respectively (with the 
latter standing for the speed of sound if thermal pressure 
dominates), and defining and Q g following eq. (7) sep- 
arately for each component, the effective Q relevant for 
the instability of the combined system is approximately 



1 



2Q- 1 



1 



(50) 



where cr g * = er g /er* and q is the dimensionlcss wave num- 
ber q = ka^/n. The first term has to be slightly modi- 
fied to take into account the dissipationless nature of the 
stars 6 , but this correction makes only a small difference 
(Rafikov 2001). The system is unstable once Q < 1. The 
most unstable wavelength corresponds to the q that min- 
imizes Q; it lies between q = 1 for cr g * = 1 and q ~ cr^ 1 
for (7 g » <C 1. Note that with <t* > cr g , the stellar disk 

6 for a dissipationless component, the q dependence in the first 
term should be replaced by q~ 1 [1 — exp (— q 2 )Io (<? 2 )] , where Io is 
the Bcsscl function of order 0. 



by itself may tend to be less unstable than the gas disk 
by itself, Q* > Q g , but through its contribution to the 
self-gravity that drives the instability, the stellar disk 
can help the gas component de-stabilize the disk. The 
combined system can be unstable for axisymmetric per- 
turbations even when each of the components has a Q 
value above unity. 

If <7 g * = 1 , the two components can be treated as one 
in the instability analysis. In order to illustrate the ef- 
fect of different values of siggs, consider, for example, 
the case of equal mass densities for the gas and stars, 
£ g = £*. When a g * = 1, the criterion for instability 
Q < 1 translates to Q g < 2 (or Q* < 2). If however 
the stars are "hotter", er g>t = 0.5 say, the criterion for 
instability (obtained at q ~ 1.6) becomes Q g < 1.42 (or 
Q* < 2.84), so the error made in Q c by ignoring the 
higher velocity dispersion of the stars (namely 2 versus 
1.6) is about 30%. In this case, the instability driven 
by the cooler gas component is only slightly affected by 
the "hotter" stellar component. When the stellar disk is 
much "hotter" than the cold gas, cr g » -C 1, the instability 
criterion becomes Q g < 1, so the stars become part of 
the stabilizing component of the system. 

In the solar neighborhood, (according to Binney & 
Tremaine 2008, §6.2.3), the gas fraction is ~ 25%, the 
stars are much "hotter" than the gas, a g ^ ~ 0.18, and 
the separate Q values are Q g ~ 1.5 and Q* ~ 2.7. This 
gives a combined value of Q ~ 1.2, indicating that the 
solar neighborhood is stable for axisymmetric perturba- 
tions because of the high velocity dispersion of the dom- 
inant stellar component. However, it is unstable for non- 
axisymmetric perturbations, only 20% variations in the 
surface density or the sound speed of the gas can generate 
axisymmetric instability. 

We note that disk stabilization can be helped by the 
tendency of the gas temperature in the smooth, warm 
disk component not to drop significantly below 10 4 K, 
corresponding to a speed of sound ~ 15kms _1 . This is 
because atomic cooling is ineffective at lower tempera- 
tures, and because supernova and stellar feedback tend 
to heat the gas to such temperatures (Wolfire et al. 2003). 
With V ~ 200 kms -1 , such a minimum speed of sound 
would not allow Q to be kept below unity once 5 < 0.13. 
This implies that once the cold disk becomes less than 
20% of the total baryonic mass, the disk tends to become 
stable against the formation of giant clumps. 

Could the z ~ 2 thick disks of the SFGs evolve to the 
thick disks of today's spiral galaxies? In the massive disk 
galaxies at low redshift, the thick disk is typically about 
a quarter of the total disk mass, and in smaller disks it 
could be as high as one half of the disk mass (Yoachim 
& Dalcanton 2006). If at z ~ 2 half the baryonic mass 
is in the disk (5 ~ 0.3) that becomes the thick stellar 
disk of today, and if, say, the total baryonic mass has 
doubled since then and all of it has settled in a new thin 
disk, then today's thick disk fraction is expected to be 
about a third, marginally consistent with the observed 
fractions. On the other hand, the velocity dispersion in 
the high-z stellar disks might be too high to be consis- 
tent with today's thick disks. Many of them may end up 
as SO galaxies, or as elliptical galaxies if they go through 
significant mergers. The simulations described in §12 
(Ceverino et al. 2009) indeed demonstrate that the evo- 
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Fig. 2. — Gas surface density of a galaxy at z ~ 2.3 from a high- 
resolution cosmological simulation (Ceverino et al. 2009). The box 
side is 15 kpc. The color code is log surface density in units of 
Mq pc~ 2 . Only gas with density n > 0.1cm -3 is shown. The 
face-on view (top) shows an extended disk broken into several gi- 
ant clumps and sheared perturbations, similar to many observed 
SFGs. The edge-on view demonstrates that this is a well defined 
disk, resembling the observed "chain" systems. The disk is fed by 
streams, not emphasized in this picture because they arc largely 
below the threshold density for the plot. With 5 ~ 0.26, this sys- 
tem is in a steady state at the late stages of the clump-instability 
phase. This galaxy, in a halo of 4 X IO^Mq, is smaller than the 
typical SFG, but it resembles the typical appearance of SFGs and 
it demonstrates the general behavior predicted by our theoretical 
analysis. 

lution into an SO-likc configuration is likely. However, 
it is clear that not all the massive high-redshift SFGs 
end up as SO's today, because the mean comoving num- 
ber density of SFGs ( £ 10~ 4 Mpc) is higher than that 
of today's SO's (< l(T 4 Mpc). 

12. COSMOLOGICAL SIMULATIONS 

Figure 2 provides a sneak preview of the results from 
zoom-in cosmological simulations that are analyzed in 



more detail in a companion paper (Ceverino et al. 2009). 
Shown are the first maps of gas density in a properly 
resolved galaxy at z ~ 2.3, which clearly show a disk 
broken into giant clumps and sheared features. They are 
brought here for the sake of demonstrating the validity 
of the simple analysis presented above. 

These zoom-in cosmological simulations are performed 
about individual galaxies using a Hydrodynamic Adap- 
tive Refinement Tree code (Kravtsov et al. 1997; 
Kravtsov 2003) with a spatial resolution better than 
70 pc (in physical units). The code incorporates the rele- 
vant physical processes, including atomic and molecular 
gas cooling and photoionization heating, star formation, 
metal enrichment, and feedback from stars and super- 
novae, as described in Ceverino & Klypin (2009). The 
unique feature of this code for the purpose of simulating 
disk instability is that it allows the gas to cool down to 
~ 100 K. This is a key to resolving the turbulence Jeans 
mass and permitting the disk fragmentation into giant 
clumps, not seen before in cosmological simulations. 

The galaxy shown in Fig. 2 is one of a few simu- 
lated galaxies, most of which showing similar features, 
although they are on the small side of the massive SFGs 
observed (as they were originally selected for another 
purpose, to match a virial mass of 10 12 M© at z = 1). In 
the specific case shown, at z = 2.3, the halo virial param- 
eters are M v ~ 4 x 10 n Af Q and R v ~ 70 kpc. The disk 
extends to R d ~ 6 kpc with a rather flat rotation curve of 
V ~ 180 kms -1 . The disk mass is M d ~ 1.1 x 10 10 M Q , 
in comparable fractions of gas and stars. The stellar 
spheroid is M sp h ~ 1.0 x 10 10 M Q , and the total mass 
within R d is M tot ~ 4.2 x 10 10 M Q , so the disk-to-total 
ratio is 5 ~ 0.26. This value is near the predicted steady- 
state value for 7 £ 0.33. 

The map clearly reveals a disk configuration dominated 
by a few giant clumps, containing a fraction a ~ 0.12 
of the disk mass. By inspecting the time evolution of 
this galaxy, we deduce that most of the clumps seen 
in the map were formed in situ in the disk, and only 
two of them came from outside. These two are embed- 
ded in dark-matter halos while all the other clumps are 
"naked" baryonic clumps. The edge-on view shows that 
the clump deviation from the mid-plane is comparable 
to and smaller than the intrinsic disk thickness. This is 
both evidence for in situ clump formation and a typical 
characteristic of observed "chain" galaxies (Elmegreen & 
Elmcgrecn 2006). With the bulge comparable to the disk, 
with the stellar component already quite developed, and 
with a non-negligible dumpiness of 7 ~ 0.2 in the incom- 
ing streams, this galaxy seems to be at a late stage of its 
giant-clump phase. The value of Q for the combined gas 
and stellar disk is below unity in parts of a ring of ra- 
dius ~ 5 kpc, but it is above unity for the disk as a whole. 
This may indicate that this system is at a stage where it is 
gradually becoming stable against axisymmetric pertur- 
bations by its old stellar component. Indeed, by the end 
of this simulation at z = 1.3, the galaxy becomes dom- 
inated by its stellar spheroid with only a minor, rather 
stable gas disk. Note, however, that the comparison of 
the linear instability analysis with the simulated clumpy 
disk should be interpreted with care, taking into account 
the fact that the linear analysis is valid at the onset of in- 
stability while the quantities quoted above are measured 
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from the highly perturbed disk. 

Convergence tests that are described in detail in Cev- 
erino et al. (2009) demonstrate that the in situ clumps 
are a result of a real physical fragmentation process and 
not a numerical artifact. In order to prevent artificial 
fragmentation due to an unresolved Jeans scale, we im- 
plemented in the code a standard procedure that im- 
poses a pressure floor (not by means of a temperature 
floor), which prevents the Jeans length from falling be- 
low Aj = A^jA, where A is the cell size of the numerical 
grid. The standard value is Nj = 4 for simulations that 
last a dynamical time (Truelove et al. 1997). To verify 
the minimum value of Nj that prevents artificial frag- 
mentation in a disk that evolves for several dynamical 
times, we repeated our simulation for 400 Myr with the 
same resolution but using ATj values that range from 1 to 
10. We found that the mass in clumps and the number of 
clumps converge to constant values for ATj = 7—10, while 
there is an overproduction of clumps that worsens with 
decreasing JVj in the range iVj = 1 — 7. This indicates 
that the clumps represent a real physical phenomenon 
once Nj > 7, so we adopt Nj = 7 for our default analy- 
sis. 

? have reported independently and immediately after 
us a similarly clumpy disk in their own simulation of a 
galaxy in a cosmological background, using an AMR code 
very similar to ours and with a comparable resolution. 
Applying A^j = 4 for the pressure floor, they indeed find 
more and typically smaller clumps, as expected from our 
numerical tests. We conclude that while their simulation 
may be partly subject to numerical fragmentation, they 
basically detect the same real physical phenomenon that 
we see in our simulations. 

13. SUMMARY AND DISCUSSION 

We studied the Toomre instability of high-redshift, gas- 
dominated, massive, thick, galactic disks as they are in- 
tensely fed by cosmological cold streams. Defining the 
disk fraction 5 = Md/M tot (Rd), an unstable disk self- 
regulates itself to Q ~ Q c ~ 1 with a velocity dispersion 
a r /V ~ Q c 5/V3. The disk develops transient elongated 
sheared features, and fragments into a few in situ bound 
massive clumps, each a few percents of the disk mass 
and together involving a fraction a ~ 0.2 of the disk 
mass. The turbulence is largely maintained by the inter- 
nal gravitational interactions within the perturbed disk. 
In particular, the encounters between the disk clumps stir 
up velocity dispersion on a timescale t cnc ~ 2.lQ 2 a~ 1 td, 
where td = — 50 Myr, which is sufficient by itself 
for matching the natural timescale for turbulence decay, 
idis ^ lAQ~H d , once Q ~ 0.67 and a ~ 0.2. 

The same gravitational encounters and dynamical fric- 
tion make the giant clumps migrate to the center on a 
timescale t m i g ~ 2.1Q 2 S~ 2 td and grow a bulge. The as- 
sociated evacuation timescale for the entire disk mass 
is iovac — 10. 5a~^ 2 1 Q 2 5~ 2 td, which is comparable to the 
timescale for mass inflow due to torques of the transient 
features. 

The cosmological streams feed baryons to the galaxy 
on a timescale t acc ~ 44t id; with r ~ 1 at z = 2 
(and varying from 0.4 to 2.5 between z ~ 1 and 9). 
The smooth component of the incoming streams, includ- 
ing small clumps, replenishes the evacuating disk, while 
the external massive clumps associated with the streams 



merge to the bulge. If the galaxy is fed by streams that 
contain less massive clumps than average, 7 < 0.33, the 
system settles into a near steady state with 8 ~ 0.3, 
where the input by streams and the transport from disk 
to bulge maintain a constant bulgc-to-disk ratio near 
unity. 

In galaxies where the incoming streams are more 
clumpy than critical, 7 > 0.33, the external massive 
clumps merge into a dominant spheroid, 8 < 0.3, and the 
dense clumps stir-up the turbulence in the disk to levels 
that stabilize the disk. The dependence of the instability 
on the degree of dumpiness in the streams introduces a 
bimodality in the galaxy properties starting already at 
z > 3. Streams smoother than average lead to extended, 
unstable disks with giant in situ clumps that form stars 
at a high rate, while more clumpy streams help building 
compact, massive bulges and stir up the disks to stabi- 
lization followed by suppressed SFR, as follows. 

In about half the halos of ~ 10 12 M© at z > 1, where 
the streams are relatively smooth, the rapid accretion at 
high z leads to a dense, gas-dominated disk with 8 £ (3 ~ 
0.6. Such a disk is wildly unstable; it grows sheared per- 
turbations and fragments into a few giant clumps, each 
comprising a few percent of the disk mass. The clump 
interactions self-regulate the disk in a marginally unsta- 
ble state, with Q ~ Q c . The clumps migrate inward in a 
couple of rotation times and form a spheroid. The com- 
bined effects of penetrating streams, disk evacuation and 
spheroid growth make 8 approach a near steady state 
within one Hubble time. They maintain the giant-clump 
phase with 8 ~ 0.3 and a r /V ~ 0.15 — 0.3 for several 
Gigayears, till well after z ~ 2. During this phase, as the 
accretion rate gradually slows down, the bulge-to-disk ra- 
tio grows slowly, the velocity dispersion oy/V declines ac- 
cordingly to keep Q ~ Q c , the clumps get slightly smaller 
in proportion to S 2 , and the clump migration and disk 
evacuation slow down by a similar small factor. 

In the other half of the massive galaxies, where the 
clumpy component of the streams is higher than aver- 
age, the dense clumps stir up turbulence in the disk at 
a level that can stabilize the disk against the formation 
of giant clumps once the bulge is more massive than the 
disk. Even if the system starts as an unstable disk form- 
ing stars at a very high redshift, once fed by clumpy 
streams it settles within a Hubble time to a phase of 
growth where the disk is stabilized. In this case, most 
of the star formation and the most rapid growth of the 
spheroid occurred during the first Gigayear or two. Once 
the disk is stable, it does not form new clumps, and when 
the remaining disk clumps have consumed most of their 
gas or disappeared, the SFR is substantially suppressed. 7 

The proposed picture has interesting implications on 
the sizes of galaxies at high rcdshift. After the system 
have grown a bulge comparable to the disk, the disk is 
expected to be twice as extended as implied by the stan- 
dard spin-parameter argument, because it acquired the 
angular momentum lost by the material that migrated to 
the bulge. Another reason for the extended appearance 
of the disk is that the coherent streams tend to join the 
disks in their outer parts. This makes the generation of 

7 This is yet another context where clumpy streams help the 
quenching in massive halos where the potential well is deep (Dekel 
& Birnboim 2008; Khochfar & Silk 2009). 
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clumps and the resulting star formation most efficient in 
an outer ring. In turn, the typical high-redshift spheroids 
are expected to be compact because the timescales for in 
situ clump migration, as well as the timescale for incom- 
ing mergers, are comparable to the star formation time in 
the clumps, making the coalescence into the bulge highly 
dissipative. 

The disk giant clumps are the cites of intense star for- 
mation, in an overall rate that follows the gas accretion 
rate of ~ lOOMQyr" 1 at z ~ 2. The apparent star- 
formation efficiency in the clumps is 77 ~ 0.1 compared to 
star formation on a dynamical time. For a local efficiency 
at a standard level of 1%, as observed in star-forming 
molecular clouds, the star formation has to be confined 
to a cuspy core or subclumps that are denser than their 
host giant clumps by one to two orders of magnitudes. 

The giant clumps are not expected to be disrupted by 
supernova feedback, but a significant fraction of their gas 
may be expelled back to the disk by radiative stellar feed- 
back (Murray et al. 2009). This would not change the 
predicted steady state in a qualitative way, but it would 
slow down the bulge growth, leaving mergers as the dom- 
inant mechanism for spheroid formation. The distribu- 
tion of stellar ages within the individual disk clumps can 
provide an observational constraint on the actual level of 
clump disruption prior to their coalescence with the cen- 
tral bulge. If the clumps indeed underwent an effective 
gas removal, the spread of ages within each clump is not 
expected to exceed ~ 100 Myr, and the outflows from the 
clumps should be detectable observationally. 

If clump disruption is not too effective during their 
migration, the efficient formation of bulges by clump mi- 
gration in the high-redshift disks helps reconciling the 
models of galaxy formation with the observed abundance 
of stellar spheroids. The cosmological major-merger rate 
seems insufficient for the purpose (e.g., Jogce 2008; Dckcl 
et al. 2009; Bundy et al. 2009). Indeed, not having 
enough mergers in the simulations, most semi-analytic 
models of galaxy formation had to assume a high rate of 
bulge growth from disk instabilities in order to match the 
observed abundance of spheroids (Cattaneo et al. 2008; 
Parry et al. 2009). Our analysis spells out the origin of 
this instability. 

As proposed by Elmegreen et al. (2008b), the massive 
black holes observed in spheroids already at high redshift 
can naturally originate from seed black holes that have 
formed by massive-star coalescence in the dense stellar 
clusters at the centers of the giant clumps. These black 
holes were shown to migrate with the clumps into the 
central spheroid, and to reproduce the observed black 
hole to spheroid mass ratio of ~ 10~ 3 . The AGN feed- 
back associated with the central black hole can add yet 
another quenching mechanism (Cattaneo et al. 2009), 
which also becomes more effective as the bulge grows. 

If the properties of the streams feeding a given galaxy 
vary in time, and in particular if the degree of clumpincss 
in them evolves, the galaxy may go through transitions 
from an unstable disk-dominated configuration to a sta- 
ble bulge-dominated state and vice versa. However, after 
z ~ 2, the recovery from a bulge-dominated system back 
to an unstable disk takes several Gigayear and may never 
materialize. 

A systematic change in the stream properties is ex- 
pected after z ~ 1, where the cosmological accretion rate 



becomes slower and the smooth cold streams no longer 
penetrate very effectively through the shock-heated me- 
dia in massive halos of ~ 10 12 M Q or higher (Keres et al. 
2005; Dekel h Birnboim 2006; Cattaneo et al. 2006; 
Ocvirk et al. 2008). Once the accretion cannot replenish 
the disks on a time scale comparable to the timescale 
for the disks to turn into clumps and stars, the galaxies 
become star dominated and eventually stable against ax- 
isymmetric modes. Supernova and stellar feedback can 
then add to the stabilization of the gas disk. The com- 
mon late disks form predominantly in halos below the 
threshold mass of ~ 10 M Q , and not necessarily by 
narrow streams (Birnboim & Dekel 2003; Binney 2004; 
Keres et al. 2005; Dekel & Birnboim 2006; Birnboim et al. 
2007). They evolve secularly through non-axisymmetric 
modes of instability associated with quiescent star for- 
mation. The turbulent high-redshift disks may end up 
in today's thick disks and SO galaxies, or in ellipticals 
through mergers. The spheroids continue to grow ac- 
cording to the standard scenario; bulges develop by slow 
secular evolution in the disks (Bureau & Athanassoula 
2005; Athanassoula 2008), and all spheroids, including 
today's giant ellipticals, grow by minor and major merg- 
ers. 

We conclude that the typical high-redshift massive 
galaxies are in a phase of evolution that does not have a 
common parallel in low-rcdshift galaxies. On one hand, 
the intense and deeply penetrating coherent, cold gas 
streams keep the disk gas rich and thus drive a wild in- 
stability with giant clumps forming stars at a high rate, 
and maintain this phase in steady state for cosmolog- 
ical times. In parallel, a high degree of dumpiness in 
some of the high-flux streams could stabilize the disks 
against the formation of giant clumps and form massive 
spheroids with low SFR already at high redshift. We 
note, however, that there are rare cases at low redshift of 
relatively gas-rich galaxies that somewhat resemble the 
perturbed clumpy appearance of the high-redshift SFGs 
(e.g., NGC 4303, Boissier et al. 2003). They may be 
scaled down versions of the SFGs, possibly unstable to 
axisymmetric modes but with much lower gas densities 
and SFR, and probably fed by overly intense accretion 
compared to the average at low redshift. 

The quenching of disk star formation by a domi- 
nant stellar bulge, which can be termed "morphological 
quenching" , is addressed and demonstrated using simu- 
lations in Martig et al. (2009). It can explain the ex- 
istence of red-and-dead early-type galaxies in the field, 
that is in halos below the critical mass for virial shock 
heating, above which the quenching can be explained in 
other ways involving termination of gas supply(Dckcl & 
Birnboim 2006). Morphological quenching predicts the 
presence of non-negligible stable gas disks in some of the 
field early-type galaxies. 

Preliminary tests indicate that the theoretical frame- 
work proposed here is in general agreement with the evo- 
lution of galaxies in hydrodynamical cosmological simu- 
lations of appropriate resolution (Ceverino et al. 2009), 
as well as with the observations of high-redshift galaxies. 
However, the comparison of the linear instability analy- 
sis with the evolved, nonlinearly perturbed systems, as 
simulated or as observed, should be performed with care. 
Our current analysis is meant to provide a simple, basic, 
theoretical framework, but the exact numerical values 
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quoted for the parameters, averaged over the whole disk, 
should not be interpreted too literally. The simulations 
allow us to predict the observable appearance of the cold 
streams that drive galaxy formation at high redshift, ei- 
ther in emission as Lyman-alpha Blobs (Dijkstra & Loeb 
2009; Goerdt et al. 2009), or in absorption as Lyman- 
limit systems or Damped Lyman-alpha systems (Dekel 
et al. 2009). 
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APPENDIX 

STEADY STATE: APPROXIMATE SOLUTIONS 

In the limit 4& 3 /(27c 2 ) <C 1, the expression in eq. (38) becomes u ~ c 1 / 3 . In our fiducial case at z ~ 2, this limit 
corresponds to (1 — j) 2 ^> 4/27, namely smooth streams, leading to steady state with the disk more massive than the 
bulge. For < 7 < 0.33 (where 0.33 > 5 SS > 0.26), this approximation underestimates S ss by a few percent. 

In the other limit, c 2 /6 3 <C 1, there is an approximate solution of the sort <5 SS = (c/b)(l — e) with e < 1; to first order 
in e it is e ~ c 2 6~ 3 /(l + 3c 2 6~ 3 ). In our fiducial case at z ~ 2, this limit corresponds to (1 — j) 2 <C 1, namely clumpy 
streams, leading to steady state with a dominant bulge. This approximation overestimates 6 SS by a few percent for 
7 > 0.3 (S BB < 0.27). 

The general behavior of the solution of eq. (32) can be evaluated analytically in the fiducial case when noticing that 
the second term involving i acc becomes negligible near 6 ~ /3(1 — 7), that is ~ 0.33 for the fiducial case, not far from 
the steady-state value. With oc t, the simplified equation 5 ~ — becomes <5/<5 3 ~ — 6.33 a.aA^ t -1 , and its 
solution is 

S- 2 - Sr 2 ~ 12.7 a 2 A^ 1 ln(t/*m) . (Al) 

According to this, with <5i n ~ 0.5, the system is expected to evolve to near 5 SS by t ~ 2<j n , as seen in the exact solution, 
Fig. 1. 

If we assume instead that id is the same at all times, the equation is replaced by S/S 3 ~ —2.0 a^t^g and the solution 
becomes 

S- 2 - Sr 2 ~ 4.0 a 2 t£ (t - t in ) . (A2) 

Again, the system is expected to evolve rapidly to near steady state, crudely reproducing the behavior in Fig. 1. 

The term neglected in eq. (32) tends to add a negative contribution to 5 for S > /3(1 — 7), and a positive contribution 
for smaller values of <5. This speeds up the evolution when the system is far form the steady state value, and it may 
either speed it up or slow it down near steady state. 

COMMENTS ON THE ORIGIN OF PRESSURE SUPPORT 

Since the gas cools rapidly to 10 4 K, the thermal pressure cannot support a clump of a r ~ 30kms _1 against 
gravitational collapse and star formation on a free-fall time. The required pressure support must be due to velocity 
dispersion. Can the gravitational interactions in the perturbed disk drive the turbulence inside the clumps? We 
address this issue using simulations elsewhere, and bring only two preliminary considerations here. 

The dissipation rate of turbulence inside a rather uniform gas clump is 

E dis c 1 -^lc2AM c a 2 t~\ (Bl) 

assuming a radius 0.5i? c for the clump, and using in the second equality R c /a c ~ t d with t<j the disk dynamical time. 
The associated timescale is tdis — 0.6-Rc/c c , i.e., comparable to the dynamical time of the clump, t c ~ 0.5id- In 
comparison, the rate of gravitational work done on the clump as it migrates a radial distance i?d is 

E mig ~ G f totMc ~ 0.5 M c V 2 Q~ 2 S 2 . (B2) 

fid t-mig 

The ratio of the two is -Edis/^mig — 1.7 Q , which is of order unity for Q ~ 1 and about a third for our fiducial 
Q = 0.67. This implies that the gravitational power associated with the migration is in principle enough for balancing 
the turbulence decay inside the clumps. The actual mechanism for pumping up the internal energy in the clumps 
could in principle be clump encounters, shear and tidal interactions with the transient perturbations, yet to be studied 
in detail. 

The decay of turbulence within the clumps may actually be slower than implied by eq. (Bl) and thus easier to 
balance. For example, the dissipation rate naturally slows down in proportion to the gas fraction as the gas turns into 
stars. The fragmentation to dense subclumps, advocated in §8, may have a similar effect by itself. If the giant clump 
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fragments to -/V su b subclumps of equal mass in which the gas density is n su b, the turbulence decay rate would change 
in proportion to the total cross section for subclump collisions, 

E dis <x iV sub ^ rx (-M . (B3) 

If ?^sub/^c ~ 36, to allow rj ~ 0.01 in the subclumps (§8), we obtain that the turbulence decay would slow down as long 
as iV c < 10 3 . This implies for 10 9 M Q clumps that the dissipation timescale would be longer than the clump dynamical 
time as long as the actual gas clouds forming stars are bigger than 1O 6 M0. 
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